Evolution of the Fermi Surface across a Magnetic Order-Disorder Transition in the 
Two-Dimensional Kondo Lattice Model: A Dynamical Cluster Approach 
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We use the dynamical cluster approximation, with a quantum Monte Carlo cluster solver on 
clusters of up to 16 orbit als, to investigate the evolution of the Fermi surface across the magnetic 
order- disorder transition in the two-dimensional doped Kondo lattice model. In the paramagnetic 
phase we observe the generic hybridized heavy fermion band structure with large Luttinger volume. 
In the antiferromagnetic phase, the heavy fermion band drops below the Fermi surface giving way to 
hole pockets centered around k = (7r/2,7r/2) and equivalent points. In this phase Kondo screening 
does not break down but the topology of the resulting Fermi surface is that of a spin-density wave 
approximation in which the localized spins are frozen. 
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In its simplest form, the Kondo lattice model (KLM) 
describes a lattice of spin 1/2 magnetic moments cou- 
pled antiferromagnetically via an exchange coupling J 
to a single band of conduction electrons and is believed 
to capture the physics of heavy fermion materials such as 
CeCue- The huge mass renormalization can be attributed 
to the coherent superposition of individual Kondo screen- 
ing clouds and the resulting metallic state is character- 
ized by a Fermi surface (FS) with Luttinger volume con- 
taining both conduction and localized electrons. At con- 
stant density, the J dependence of the inverse effective 
mass — or the coherence temperature — has been argued 
to track the single ion Kondo scale [H,!!]. Competing with 
Kondo screening, the localized spins interact indirectly 
via magnetic polarization of the conduction band. This 
RKKY scale dominates at low values of the exchange cou- 
pling and is the driving force for the observed magnetic 
order-disorder quantum phase transitions in those mate- 
rials . The nature of this phase transition is of current 
interest following experimental results suggesting a sud- 
den change in the FS topology at the quantum critical 
point (QCP) for the heavy fermion metal YbRh2Si2 0- 
Driving this system from the non-magnetic heavy fermion 
metallic phase to the antiferromagnetic (AF) metallic 
phase causes a rapid change in the low temperature Hall 
coefficient which is extrapolated to a sudden jump at 
T = 0. Since the low-temperature Hall coefficient is re- 
lated to the FS topology the results are interpreted as 
showing a sudden reordering of the FS at the QCP from 
a large FS, where the local moment impurity spins are 
included in the Luttinger volume, to a small FS where 
the impurity spins drop out of the FS volume. This sce- 
nario lies at odds with the Hertz-Millis description of the 
quantum phase transition [sl, 0] and has triggered alter- 
native descriptions 0, S] • In this study it is exactly this 
issue which we investigate using the Kondo lattice as our 
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FIG. 1: (a) Definition of the real space basis vectors, unit 
cell and orbital indices (b) The reduced Brillouin zone and 
k-space patching for a cluster with N^^ — 4 



model system: 



F = ^e(k)4^,Ck.. + J^Sf.Sf 
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with c\. ^ creating a conduction electron on an extended 
orbital with wave vector k and a z-component of spin 
(J =t,i. The spin 1/2 degrees of freedom, coupled via 
J, are represented with the aid of the Pauli spin ma- 
trices (T by S? = \ gf c\ ^ds^s'C'i.s' or the equivalent 
definition for S/ using the localized orbital creation op- 
erators fl^. The KLM forbids charge fluctuations on the 
/-orbitals and as such the constraint of one electron per 
localized orbital must be included. To avoid particle-hole 
symmetry at half-band filling, we opt for a dispersion re- 
lation: e(k) = —2t[co^{kx) + co^{ky)] — 2t'[Q.os{kx + ky) + 
(zos{kx — ky)] where t and t' are the ffist and second near- 
est neighbor hoping matrix elements respectively. We use 
t' /t = —0.3 throughout our calculations. 

Method - To solve the above model, we apply the 
dynamical cluster approximation [13] (DCA) with a 
Hirsch-Fye quantum Monte Carlo (QMC) algorithm as 
cluster solver. The DCA relies on a /c-space coarse 
graining. The Brillouin zone is divided into Nc patches 
with center K and the original fc- vectors are given by: 
k = K -\- k. The DCA boils down to satisfying mo men- 
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turn conservation only between patches. As can be seen 
within a skeleton expansion, this approximation yields 
a self-energy which is a functional of the coarsed-grained 
Green function: G{K,iUn) = ^Y^kepatch^^^^^^^^ri)- 

With G{k,iuOn) = {G^^ik.iujn) - T.[G{K,iuOn)\)~^ we 
obtain a self-consistent equation for G. Each iteration 
step requires the solution of Nc impurities embedded in 
a fermionic bath given by: Q^^ = G~^ + S. 

Our particular implementation of the DCA allows for 
spontaneous antiferromagnetic symmetry breaking which 
enables us to not only localize the magnetic transition but 
also to carry out simulations within the magnetically or- 
dered phase. The setup of our calculation is shown in 
Fig. [1^. Broken lattice symmetries are taken into consid- 
eration by defining a unit cell of two lattice sites contain- 
ing, in total, two conduction electron orbitals and two 
localized /-electron orbitals. Within this unit cell, the 
up and down bath Green functions are allowed to take 
different values thus generating spin symmetry broken 
solutions. This is the smallest cluster, which we denote 
by N^^ = 1, with which we can capture AF ordering. 
Translational symmetry is now only assumed for the new 
basis vectors ai and a2 giving the magnetic Brillouin zone 
(MBZ) in momentum space. DCA patching of the MBZ 
is demonstrated in Fig. [TJd for a cluster of size N^^ = 4. 
Importantly, the DCA self-consistent equation for the lat- 
tice Green function becomes a 4 x 4-matrix equation re- 
flecting the four orbitals of the unit cell. In this study we 
use both N^^ = 1 and N^^ = 4 clusters. Although, for 
instance, with A^^^ = 1 the coarse-grained self-energy 
has no momentum dependence a full fc-dependency of 
the lattice Green function enters via the non-interacting 
single-particle Green function Go{k, iujn), thereby allow- 
ing the calculation of a detailed FS regardless of the clus- 
ter size. Our implementation of the Hirsch-Fye impurity 
algorithm [ijj to solve the cluster problem follows pre- 
cisely the ideas introduced in Ref. [12]. 

Results - The DCA is justified if the underlying physics 
is driven by the frequency dependence of the self-energy 
as opposed to its momentum dependence. To check this 
we have carried out simulations at the particle-hole sym- 
metric point {t' /t = and (ric) = 1) and compared the 
results with the lattice QMC simulations of Ref. [l2[ • Al- 
though the DCA at A^^^ = 1 grossly overestimates the 
coupling at which the magnetic order-disorder transition 
occurs, it reproduces the single particle spectral function 
both in the paramagnetic (PM) phase and, most impor- 
tantly, in the AF phase remarkably well. This is demon- 
strated in Fig. [2] for the AF phase. Hence, the DCA 
is able to capture the delicate interplay between Kondo 
screening and magnetic ordering [12]. The sign problem 
which restricts finite sized lattice QMC studies of the 
KLM to the particle-hole symmetric case is not severe 
in the DCA approach close to half-filling and on cluster 
sizes up to 16 orbitals. Fig. [3] maps out the ground state 
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FIG. 2: (color online) Comparison of the single particle spec- 
tral function of the KLM at half- filling (t'/t = and {uc) = 1) 
and with coupling J/t = 1.2 from (left) 12 x 12- lattice QMC 
results as obtained with the projective auxiliary field algo- 
rithm of Ref. [m and (right) DCA results with cluster size 
N^^ = 1 and /3t = 40.0 
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FIG. 3: (color online) Bottom: Magnetic phase diagram of the 
hole-doped KLM showing simulation results for the staggered 
magnetization m{ (color-coded circles) as a function of cou- 
pling J/t and conduction electron occupancy (ric). Shading 
of AF and PM regions is a guide to the eye. Here t' /t = —0.3 
and the calculations are carried out with the N^^ = 1 cluster. 



Top: 



as a function of (ric) at constant coupling J/t — 1 



magnetic phase diagram of the KLM as a function of cou- 
pling J/t and conduction band hole-doping. Here, we are 
interested in ground state properties, so the choice for the 
inverse temperature (3t must be large enough to ensure 
we are below the smallest scale in the problem: the coher- 
ence scale and/or the RKKY scale. Since the coherence 
scale decays exponentially with J / 1 we are limited to val- 
ues of J/t > 0.8. This restriction arises since the compu- 
tational time required by the QMC cluster solver scales as 
{(3N^^Y [ill- general grounds, we expect the onset 
of magnetism at small values of J/t since in this region 
the RKKY scale set by J^xIq^ ^ = 0) dominates over the 
Kondo scale, Tk ~ e~^l'^ . Here, xIq?^ = 0) corresponds 
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FIG. 4: (color online) DCA single particle spectrum for simu- 
lation A (see Fig.E)) in the PM region ( J/t = 1.0, (uc) = 0.855 
and pt = 40.0). The bottom right plot is a close-up of the 
main plot for energies around the Fermi-energy along the path 
(0,0) to (7r,7r). The bottom left plot shows the resultant 
Fermi- surface, plotted here via a large- A/" mean-field model- 
ing of the DCA data (see text) 



FIG. 5: (color online) DCA result for the single particle spec- 
trum with N^^ = 4, J/t = 1.0, (ric) = 0.977 and pt = 40.0 
in the AF region (point B in Fig. [3}. The bottom right plot 
is a close-up of the main plot for energies around the Fermi- 
energy along the path (0,0) to (7r,7r). The bottom left plot 
shows the resultant Fermi-surface plotted here by fitting to 
the DCA spectrum using the mean-field form of Eq. [2] 



to the spin susceptibility of the conduction electrons. It 
is interesting to note that within the DCA approxima- 
tion on cluster sizes up to N^^ = 4 the value of J/t 
at which antiferromagnetism sets in at half-band filling, 
ric = 1, is only slightly affected by the choice of t' /t. On 
the other hand, the charge gap, which is known to scale 
with J/t in the particle- hole symmetric case {f /t = 0) 
[l2| , is very much suppressed away from the particle- hole 
symmetric point. In particular, at t'/t = —0.3 our re- 
sults are consistent with an exponential scaling of the 
charge gap [13]. Upon doping we observe a magnetic 
metallic state which as a function of decreasing coupling 
J/t progressively dominates the phase diagram. Finally, 
within our numerical accuracy, the staggered magneti- 
zation vanishes smoothly at the magnetic order-disorder 
transition thus lending support to a continuous transition 
(see Fig. [31 top panel). 

Having established the magnetic phase diagram, we 
turn to our primary concern, namely the evolution of 
the FS across the magnetic order-disorder transition. 
Figs. [4] and [5] plot the single particle spectral function, 
A{k,u;) = -ImG^f (fc,cj) in the PM and AF metallic 
phases corresponding to the points A and B respectively 
in the phase diagram (Fig. [3]). Analytic continuation to 
real frequencies was achieved with a stochastic analytic 
continuation method [3, [lH . 



The DCA spectral function in the PM phase (Fig. [4]) 
is characterized by an extremely fiat heavy fermion band 
crossing the FS in the vicinity of k = (tt, tt) and equiva- 
lent points. In the AF phase (Fig. [5|), again at J/t = 1.0 
but with N^^ = 4 and in close proximity to half-band 
filling the heavy fermion band centered around (tt, tt) 
is still present but has dropped below the Fermi en- 
ergy and hole pockets centered around the wave vector 
k = (7r/2,7r/2) arise. This result is further supported by 
calculations at J/t = 0.8, pt = 80, closer to the mag- 
netic transition {uc = 0.92) but on smaller cluster sizes 
{N^^ = 1). The continued presence of heavy bands 
around (7r,7r), although shifted to lower energies, indi- 
cates that partial Kondo screening still occurs in the 
magnetic phase. The main result of our spectral anal- 
ysis however is the topological change of the FS when 
comparing the PM and AF phases. 

Interpretation - As an aid to interpreting our DCA 
results we draw on the mean-field approximation of 
Zhang and Yu [16] which introduces the order parame- 
ters (/i^^/i,T - flifid) = mfe'^\ (4,T^i.T - 4,i^''d) = 
—rrice'^^^, with Q = (7r,7r), and (//^Ci^i + Cj = 

+ Cj t'^i't) ~ ^/ are the stag- 

gered magnetizations of the impurity spins and conduc- 
tion electrons respectively and ^ is a hybridization order 
parameter which mimics the screening of the impurity 
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spins. In this mean-field approach, which we stress is 
used here purely as a tool to aid interpretation of our 
DCA results, the constraint of unit occupation on the 
localized /-electron orbital is taken into account only on 
average by means of a Lagrange multiplier A, as opposed 
to the exact single occupancy enforced within the DCA- 
QMC cluster solver. 

The mean-field Hamiltonian, given by 
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results in a four-band energy dispersion relation, £^n(k), 
in the MBZ. 

The low temperature features of the PM phase are well 
understood in this mean-field framework by setting m/ = 
mc = in Eq. ([2]) which recovers the generic hybridized 
band structure of the large-TV approximation, £^±(k) = 

'^-^^^ ± \ (^(ek - /i - A) V ( JV)2) ' with heavy bands 
crossing the Fermi energy in the vicinity of k = (tt, tt) and 
equivalent points in excellent agreement with the DCA 
result. In this mean-field modelling of the PM phase the 
Luttinger volume is given by IVbz = (^/ + ^c)/2, 
counting both conduction band electrons and impurity 
spins. Since the mean-field model accounts very well for 
the DCA spectral function in the PM phase, with the 
exception of the spectral weight around k = (0,0) and 
cj/t -1/4 [18], we attribute the DCA result with this 
same large Luttinger volume. 

In the AF phase the DCA spectrum can again be well 
accounted for with the mean-field Hamiltonian of Eq. ([2]) 
but now setting non- vanishing staggered magnetizations, 
equal to the QMC measured observables, and using a 
non-zero hybridization value of F = 0.3 as fit parameter. 
In this fit, one band drops completely below the Fermi 
energy which itself is crossed only by the second band 
£^2(k) in the vicinity of k = (7r/2,7r/2). In this case 
the Luttinger volume is given by jVuBZ = and 
the topology of the FS is that of a spin-density wave 
approximation = 0, m/ = l,mc 7^ in Eq. ([2])) where 
the /-electrons are frozen and do not participate in the 
Luttinger sum rule. Hence, we coin the FS, obtained 
from the DCA spectrum, as small. 

We note that similar results have recently been ob- 
served in a variational Monte-Carlo approach to the KLM 
A major difference however is that the magnetic 
phase transition separating topologically different Fermi 
surfaces is of first order in the variational approach fv^ . 
In contrast, the DCA calculation presented here supports 
a continuous transition in accordance with experimental 
findings. 



Conclusions - We have presented large scale DCA cal- 
culations of the Kondo lattice model, and mapped out 
the magnetic phase diagram in the J/t versus doping 
plane. We have tested the approximation and seen re- 
markable agreement with lattice QMC methods with re- 
spect to the single particle spectral function. In partic- 
ular, the approximation captures the delicate interplay 
between Kondo screening and magnetic ordering. Across 
the magnetic phase transition, the data shows a change 
in the topology of the FS from a large FS in the PM 
state to a small FS with hole-pockets centered around 
k = (7r/2,7r/2) in the AF phase. This change in topol- 
ogy of the FS is not linked to the breakdown of Kondo 
screening, and is consistent with a jump in the Hall co- 
efficient as observed experimentally This jump in 
the Hall coefficient is not linked to a first order transi- 
tion since within our numerical accuracy the staggered 
magnetization is a continuous function of the control pa- 
rameter driving the quantum phase transition. 
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